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Abstract 

We study coarsening dynamics in the ferromagnetic random bond 
Ising model in d = 1,2. We focus on the validity of super-universality 
and the scaling properties of the response functions. In the d = 1 
case, we obtain a complete understanding of the evolution, from pre- 
asymptotic to asymptotic behavior. The corresponding response func- 
tion shows a clear violation of super-universality. Further, our results 
for d = 1,2 settle the controversy regarding the decay exponent which 
characterizes the response function. 



A system undergoes phase ordering when it is suddenly quenched from 
a high temperature to below the critical point (T < T^) at time t = 0. The 
basic feature of this process is relaxation via domain coarsening [H |2]- In 
a system without disorder, this is a dissipative and scale-free phenomenon. 
The only relevant length scale is the typical domain size, which grows as a 
power law L{t) ~ t^^^, where z is the dynamical exponent. The evolving 
system is characterized by spatio-temporal correlation functions of the order 
parameter, which exhibit scaling properties. For example, the generic two- 
time observable 0(t, s) is expected to scale as 0(t, s) = L°'{s)fo [L(t)/L{s)], 
where {t, s) are a pair of times after the quench. We would like to obtain a 
good understanding of exponents like a, and the scaling functions fo{x) PQ[2]- 

There has also been much interest in the ordering of systems with quenched 
disorder but without frustration, where the pattern of ergodicity and symme- 
try breaking below the critical point is the same as that of the pure system. 
This includes ferromagnets subjected to random external fields or with ran- 
dom exchange couplings, modeled by the random field Ising model (RFIM) 
or the random bond Ising model (RBIM) [3]. The primary effect of disorder 
is to create energy barriers slowing down the coarsening process. The depen- 
dence of these barriers on L determines the nature of the asymptotic growth 
law, which is usually logarithmic or power-law with a disorder-dependent 
exponent. 

Apart from the growth law, a key question is how disorder affects the 
scaling function foi^), namely whether it enters only through the growth 
law or if there is also an explicit dependence. In the first case, fo{x) is 
the same for the pure and disordered systems - this is referred to as super- 
universality (SU) In the second case, the scaling relation generalizes 
to 

0{t, s) = L^{s)fo [m/L{s), L{s)/L*] , (1) 

where L* is a disorder-dependent scale. Typically, a scaling function of the 
above form implies a crossover from pre-asymptotic [L{s) <^ L*] to asymp- 
totic [L(s) > L*] behavior. 

The validity of the SU hypothesis has been mainly investigated for the 
equal-time correlation function or structure factor, and it is found to hold 
in a wide variety of cases [5l |6]. Recently, the investigation of SU has been 
extended by Henkel and Pleimling (HP) [71 [8] to the autocorrelation and 
autoresponse functions in the d = 2 RBIM, again confirming its validity. So 
far, no crossover of the type implied by Eq. ([1]) has been reported, except for 
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the = 1 RFIM P]. 

In this letter, we present a comprehensive study of the response functions 
in the ordering dynamics of the d = 1 and d = 2 RBIM. In the d = 1 case, we 
obtain a complete theoretical picture of the early- and late-stage dynamics. 
We demonstrate the existence of a crossover in the scaling functions, and the 
consequent violation of SU. Furthermore, from our study of the zero-field- 
cooled susceptibility x(t, s), we find that the corresponding decay exponent 
a = 0, as in the pure system. In the d = 2 case, it is more difficult to 
make a clear statement on the issue of SU. However, the data for s) 
conclusively shows that the scaling exponent a depends on dimensionality 
d, and is consistent with the phenomenological formula a = {d — l)/{2z). 
In fact, the presence of disorder enables us to fix this relation with higher 
precision than has been previously possible in the pure case. 

Let us first present results for the d = 1 RBIM, which is defined by the 
Hamiltonian Ti = — JiOiai+i, where the spins (jj = ±1. The ferromag- 
netic couplings Ji are independent random variables, uniformly distributed 
in the interval (1 — e, 1 + e) with < e < 1. This system undergoes ordering 
only at T = 0. However, the T = quench must be avoided to prevent the 
dynamics from getting stuck in metastable states. The way out is to quench 
to a temperature high enough to drive activated coarsening, but low enough 
to inhibit the nucleation of equilibrium fiuctuations. Then, as long as L{t) is 
smaller than the equilibrium correlation length ^(T), one observes the same 
coarsening behavior as in the T = quench. 

We have performed simulations of the Glauber-RBIM with different val- 
ues of e on a system of linear size N = 10^, up to tmax = 50000 Monte Carlo 
steps (MCS). At t = 0, the system is quenched from T = oo to T = 0.05. For 
each value of e, averages have been taken over 2000 independent realizations 
of disorder and initial conditions. We have studied the growth law [Fig.[T]^a)], 
measuring L{t) as the inverse density of defects. This is obtained by dividing 
the number of sites with at least one oppositely-aligned neighbor by the total 
number of sites [TOj. 

Figure [U^a) shows the existence of three regimes. At early times, the 
growth is very rapid and independent of e. This is the time regime where all 
micro-domains seeded by the random initial condition are eliminated with 
a few fiips. Then, the interfaces get trapped in the local energy minima, 
namely on the weakest nearby bond. This is followed by the intermediate 
regime, where growth is slowed down by activated escape over the barriers. 
This regime is absent in the pure case and lasts up to the time needed to 
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overcome the largest barrier 4e, t* ~ exp (4e/T). Finally, in the asymptotic 
regime, interfaces are effectively free and L{t) grows algebraically with z = 2, 
as in the pure case. This is clearly shown by the plot of the effective exponent 
l/zes = d{\nL) /d(\nt) vs. t/t* in Fig. [T]^b). Here, the early-time behavior 
corresponds to the initial fast drop, and the intermediate regime corresponds 
to the ensuing climb of the curves. In the asymptotic regime, the effective 
exponent becomes constant at z^s = 2 for t/t* > 1. 

Denoting the duration of the early regime by to, and integrating 1/zes 
from to onward (where it depends only on t/t*), the scaling form of the 
growth law L{t) = L(to)C(t/t* ,to/t*) is found, with 



C{t/t*,to/t*) =exp 



t/t* 

dx 



to/f XZe{i{x) 



(2) 



Hence, for values of e such that the intermediate regime is long enough to 
allow for both L{t) ^ L{t*) and l/zcf[{x) small and slowly varying, one has 
a crossover from a power-law behavior (with a disorder- dependent exponent) 
to the power law of the pure case: 

(L{to)it/toy/^'^^/^\ forto<t«r, 
^ ^ \D{e/T){t/t*y/\ for t > t*. ^ ^ 

Here, l/z{e/T) is the minimum of 1/zesix) for a given value of e/T, and 
D{e/T) = L(to)[exp /^^/j. dx{xzes)^^] is the disorder-dependent diffusion con- 
stant. The two limiting power- law behaviors are denoted by solid lines in 
Fig.IIKa). 

Our next step is to check for the validity of SU in the d = 1 RBIM. 
We have calculated the autocorrelation and autoresponse functions, which 
show scaling behavior as in Eq. ([1]) with L* = L{t*). In Fig. [21 we plot 
C(t, s) = ((jj(t)(Tj(s)) vs. t/s for different values of e and s, and a fixed value 
of the ratio q = t*/s = 1.1. This plot shows an excellent data collapse. 
According to Eq. ([T]), the collapse should occur in the plot of C{t,s) vs. 
L(t)/L{s), but it is straightforward to show from Eq. ([2]) that L(t)/L{s) 
and L{s)/L* depend on the time arguments through the ratios t/s and s/t*. 
In the inset, we plot C{t,s) vs. t/s for different values of q. In this case, 
there are different scaling curves for each q, demonstrating that the scaling 
function violates SU. Further, the curves tend towards that for the pure case 
as g — 7- 0. 
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As mentioned earlier, HP [UIH], have studied the hnear (integrated) re- 
sponse function for the d = 2 RBIM. The study of two-time quantities pro- 
vides a novel testing ground for aging and SU. The response function deter- 
mines the effect on the local magnetization at the time t, due to a small, 
constant and site-dependent external field switched on during a time inter- 
val preceding t. Different response functions correspond to different choices 
of the time interval [H] . For instance, HP study the thermoremanent mag- 
netization (TRM), which corresponds to switching the field on during (0, s) 
with s < t. Here, we study the zero-field- cooled susceptibility xit,s), corre- 
sponding to the field acting in the interval The scaling behavior of 
this quantity is an unsettled issue even for pure systems. We know that the 
scaling form xit^^) ~ ■5~°/x[-^('^)/-^('^)] obeyed in the aging regime, but 
the exponent a is still a matter of controversy |Tll [12] . Our results in this 
paper settle this issue for both pure systems and the RBIM. 

There are two different arguments for a, which originate in the context of 
pure systems. The first picture is based on the idea that the response comes 
entirely from the paramagnetic spins at the interfaces [13]. Therefore, the 
response function (per spin) should decrease at the same rate as the interface 
density L~^(t), implying 



independent of d [H]. Numerical support for Eq. (jl]) has been obtained by 
measuring the TRM [HI [15] . 

In the second picture, there is another mechanism in addition to the 
paramagnetic response, whereby domains as a whole grow so as to minimize 
the magnetic energy. This picture, which is supported by some analytical [T6| 
[T7] and numerical results [TH [19] for x(t, s), yields the phenomenological 
formula: 



The different behaviors for d < 3 and d > 3 are related to the roughening of 
the interfaces |19] . 

Despite considerable numerical effort, it has not been possible to clearly 
decide in favor of Eq. (jl]) or Eq. (|5]) for pure systems. Considering that there 
is a large discrepancy between the two expressions for d = 1,2, one may 
wonder why this is so. The root of the problem lies in the use of different 
response functions. We have argued [Hj that the TRM is affected by an 



az = 1, 



(4) 




for d < 3, 
for d>3. 



(5) 
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extended crossover, which prevents the observation of the asymptotic scahng 
behavior in the simulations [2]. On the other hand, the zero-field-cooled 
susceptibility is free of this shortcoming [TT] . 

In the light of the above discussion, let us obtain the exponent a for 
the d = 1 RBIM. The plot of x{^-,s) vs. t/s in Fig. [3] displays the same 
scaling behavior and violation of SU as C{t,s) in Fig. [2J Again, the data 
sets show excellent collapse when q = t*/s is kept constant. However, the 
scaling function depends on q and approaches the pure result as g — )■ 0. The 
collapse for many different values of s in Figure [3] shows that a = for the 
d = 1 RBIM, consistent with Eq. (jS]), and as demonstrated analytically in 
the pure system [IS]. 

Let us consider next the d = 2 RBIM. First, we focus on the domain 
growth law. Huse and Henley pO] have argued that, in this case, the bar- 
rier heights scale as AE ~ T{e)L'^. The matching of the barrier size with 
the thermal energy AE = T introduces a new length scale into the prob- 
lem, L* = (T/T)^/*^. Then, there should be a pre-asjTiiptotic regime for 
L{t) <^ L*, with the algebraic behavior L{t) ~ t^^^ of the pure system. This 
is followed by an asymptotic regime for L{t) ^ L*, where growth becomes 
logarithmic L(t) ~ [ln(t/t*)]^/'^. However, it remains controversial whether 
the crossover is from algebraic to logarithmic, or from algebraic to algebraic 
with z dependent on disorder. A number of experiments on random magnets 
pi] have suggested the latter. Further, Paul et al. fB] have obtained com- 
prehensive MC data for the d = 2 RBIM, which supports the latter scenario, 
at least on numerically accessible length-scales and time scales. However, 
Cugliandolo et al. |22l |23] have argued that this is an intermediate regime. 

We have measured L{t) from the density of defects, as in the d = 1 case, 
after quenching a 1000^ system from T = oo to T = 0.1. We average over 
different disorder realizations (from 10 to 30) for each value of e/T. In Fig. HI 
we plot the effective exponent 1/zefr vs. t, similar to Fig. [U^b). The data 
sets show the existence of an early regime, between and to, followed by a 
second regime. As in the d = 1 case, to is independent of e, but somewhat 
larger. For t > to, V-Zefr becomes approximately constant over 3 decades of 
time. On the time-scales of our simulation, there is no sign of a crossover to 
a logarithmic growth, consistent with earlier studies [211 E] • The value of the 
growth exponent is consistent with the z{e/T) = 2 + ce/T behavior proposed 
by Paul et al. [6J . Without tackling the issue of what is the "truly asymptotic 
growth law" in the d = 2 RBIM, we stress that we have at least 3 decades of 
an algebraic growth law. The corresponding exponent is e/T-dependent and 
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spans the range from z = 2 to 2; ~ 9. 

Let us turn next to the response functions in the d = 2 RBIM. The SU of 
these functions will be discussed in a later publication. Here, we focus on the 
decay exponent a. HP [8j have obtained a from the TRM for a range of z- 
values. However, their results are not consistent with either Eq. (jl]) or Eq. ([5]). 
Concerned with the violation of Eq. (jl]), HP have argued that the departure 
is due to the disorder inducing fractality in the interfaces. This would modify 
Eq. to az = d — df, with df being the fractal dimensionality of interfaces. 
In order to test this, we have examined the equal-time correlation function, 
which scales as C{r,t) = f{r/L). For fractal interfaces, the short-distance 
decay (or Porod decay) of this quantity is C{r, t) ~ 1 — b{r / LY~'^f for r/L <^ 
1. We find no difference between the short-distance behavior of C{r,t) with 
and without disorder. Therefore, disorder does not result in fractal interfaces 
and Eq. @, if correct, should also hold with disorder. We conclude that the 
HP data is unexplained due to the difficulty of accessing the asymptotic 
regime with the TRM. 

We now present our results for a in the d = 2 RBIM. For each value of 
e/T, we have computed a as the exponent producing the best data collapse 
of L°'^{s)x(t, s) vs. t/s, for different values of s. The measured products az 
are plotted in Fig. [5l together with Eq. (j4]) {az = 1) and Eq. (|5]) {az = 1/2). 
Clearly, our data is consistent with Eq. (|5]), and completely rules out Eq. @. 
We stress that both z and a vary with disorder (cf. Fig. IHfor z vs. e), but 
the product az seems to be universal. 

In summary, we have undertaken a comprehensive study of autocorre- 
lation and response functions in the ordering dynamics of the RBIM. In 
d = 1, the growth exponent shows a crossover from a pre-asymptotic disorder- 
dependent value to the asymptotic value 2; = 2, as in the pure case. The cor- 
responding autocorrelation and response functions violate SU. The scaling 
exponent of x{tj s) is a = 0, consistent with az = from Eq. ([5]). In d = 2, 
after the transient regime, we see an extended regime of power-law growth 
with z being dependent on e/T. More important, we find az ^ 1/2, again 
consistent with Eq. In general, the introduction of disorder complicates 
the domain growth problem. However, it also provides us an excellent op- 
portunity to make a clear assessment of the az-relationship. We believe that 
a further pursuit of phenomenology is of secondary importance at this stage. 
Now the focus should be on a significant theoretical advance enabling us to 
ascertain the accuracy of Eq. ([5]). 
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Figure 1: (a) Domain growth law for the d = 1 RBIM. The disorder amph- 
tude e ranges from to 0.3 in steps of 0.05 (top to bottom). The straight 
hnes depict the laws L{t) ~ t^/^; and L{t) ~ L(to)(tAo)^/^^'/^^ with e/T = 3, 
to = 11 and l/z{e/T) = 0.05 (see text), (b) Plot of l/z^s vs. t/t* for the 
data sets in (a). 
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Figure 2: Main panel: The upper curve plots C{t, s) vs. t/s with q = t*/s = 
1.1 and s = 50, 200, 800. The lower curve corresponds to the pure case. Inset: 
The same plot with three values of g = 1.1, 0.27, 0.07 (from top to bottom). 




t/s 



Figure 3: Plot of x{t,s) vs. t/s for the d = 1 RBIM with q = t*/s = 
1.1, 0.27, 0.07 (bottom to top). We present data for several values of s, rang- 
ing from 50 to 30000. 
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Figure 5: Plot of az vs. 1/z for the d = 2 RBIM. The horizontal hnes are 
drawn at az = 1 and az = 1/2. 
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